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t>- ' ABSTRACT 

i Nine quasar absorption spectra at 21-cm and UV rest- wavelengths are used to estimate 

possible variations in x = a 2 g p ft, where a is the fine structure constant, g p the proton 

, g-factor and /i = m c /m p is the electron-to-proton mass ratio. We find (Ax/x)^^ tcd — 

\Q ' (0.63 ± 0.99) x 10~ 5 over a redshift range 0.23 < z a hs < 2.35 which corresponds to 

£S) | look-back times of 2.7-10.5 billion years. A linear fit against look-back time, tied to 

■ Ax/x = at z = 0, gives a best-fit rate of change of x/x = (—0.6 ± 1.2) x 10~ 15 yr _1 . 
C> We find no evidence for strong angular variations in x across the sky. Our sample is 

, much larger than most previous samples and demonstrates that intrinsic line-of-sight 

velocity differences between the 21-cm and UV absorption redshifts, which have a 
random sign and magnitude in each absorption system, limit our precision. The data 
^ C| directly imply that the average magnitude of this difference is Ai>i os ~ 6 km s _1 . 

O Combining our Ax/x measurement with absorption- line constraints on a- variation 

yields strong limits on the variation of /i. Our most conservative estimate, obtained by 

■ assuming no variations in a or g p is simply A/i/fi = (Ax/x)™ t ^ c . If we use only the 
\ four high-redshift absorbers in our sample, we obtain A/i/fi = (0.58 ± 1.95) x 10~ D , 

d • which agrees (2<r) with recent, more direct estimates from two absorption systems 

containing molecular hydrogen, also at high redshift, and which have hinted at a 
• i— i , possible //-variation, Afi/^t = (—2.0 ±0.6) x 10~ 5 . Our method of constraining A/x//z 

■ is completely independent from the molecular hydrogen observations. If we include 

the low-redshift systems, our A/x//i result differs significantly from the high-redshift 
molecular hydrogen results. We detect a dipole variation in /i across the sky, but given 
the sparse angular distribution of quasar sight-lines we find that this model is required 
by the data at only the 88 per cent confidence level. Clearly, much larger samples of 
21-cm and molecular hydrogen absorbers are required to adequately resolve the issue 
of the variation of /i and x. 
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1 INTRODUCTION 

The existence of extra spatial dimensions remains an open 
question in theoretical physics. The answer to this question 
is closely linked to attempts at super-unification, which has 
become a modern holy grail. In turn, the existence of ex- 
tra spatial or temporal dimensions may be inferred by the 
detection of spatial or temporal variations in the values of 
coupling constants, such as the fine structure constant, a 
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(see lUzar]|2003l : [Karshenboim et alJl2005l . for reviews). This 
constant provides a measure of the strength of the electro- 
magnetic interaction and is truly fundamental because it is 
dimensionless and thus independent of any system of units. 
Other such constants are the electron-to-proton mass ratio, 
/x, the proton gyromagnetic factor g p , as well as the constant 
x defined as 

x = a 2 g p ^ , (1) 

which combines all three. Apart from the question of extra 
dimensions, any confirmed detection of a variation in the 
value of any of these 'constants' would in itself be a ground- 
breaking result. For example, in some popular models this 
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variation is related to violation of the Principle of Equiva- 
lence. 

There have been several attempts to constrain 
the variation of fundamental constants either exper- 
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the O klo nuclear reactor (Shlvakhtcr 1976: Ijjun 
| 2000l: IFlambaum fc ShurvalJ l2002t lOlive et al.1 12002 
lLamoreaux^^To^eraonn2004T) provide local, i.e. earth- 
bound and relatively recent, z ~ 0, constraints. The cosmic 
microwave background provides con s traint s at z ~ 1000 
jLandau et alJ 1200 it lAvelino etafl l200lt iMartins et"al] 
2002). Further constraints m ay be obtained from B ig Bang 
Nucleosynthesis data (see e.g. lDmitriev et ai"]|2004|) . 

However, one very promising way to detect such vari- 
ations is astronomical spectroscopy of gas clouds which 
intersect the lines of sight to distant quasars. This tech- 
nique allows one to obtain constraints essentially over all 
of the intermediate redshift range, which is crucial for 
probing variations with cosmic time. Additionally, results 
achieve very high precision, due to the advent of 10-m class 
telescopes equipped with state-of-the art high-resolution 
spectrographs, in particular Keck/HIRES in Hawaii and 
VLT/UVES in Chile. 

In this type of work, the gas clouds probed are of- 
ten damped Lyman a systems (DLAs). These systems have 
high neutral hydrogen column densities (A^hi ?J 10 20 cm -2 ) 
at precisely determined redshifts, thus providing excellent 
probes of the bulk neutral atomic hydrogen in the dis- 
tant Universe. Furthermore, DLAs contain many different 
heavy elements which absorb in the rest-frame ultraviolet 
(for atomic species) or millimetre (for molecular species), 
giving rise to a multitude of absorption lines. With the latest 
generation of optical telescopes rest-frame ultraviolet lines 
can be observed with resolutions of the order of 6 km s _1 
so that redshifts can be determined to six significant dig- 
its. This is crucial as all methods using this type of data 
aim to detect small shifts induced by a possible variation 
of constants which are of different magnitude for different 

transitions. 

Using such data, iDzuba et alJ lll999l) and IWebb et alJ 

( 1999) developed the highly sensitive many-multiplet method 
and in a series of papers applied it to rest-frame ul- 
traviolet atomic quasar absorption lines, obtained with 
Keck/HIRES. Their work provided evidence that a is 
smaller at redshifts 0.2 < z < 3.7 at the fractional level 
of Aa/a = (-0-57 ± 0.10) x 10~ 5 JWebb et al] Il999t 
Murphy et"alll200lt IWebb et aklboOll l2003t iMurphv etal 



2003, 2004 1). After thorough i nvestigation IMurphv et al 
feOOll) and lMurphv et all <20Q. " 



ruled out a large range of 
known systematic effects. However the issue remains contro- 
versial after another group found no evidence for a change in 
a using VLT/UV ES spectra (Aa/a = (-0.06±0.06) x 10~ 5 
for 0.4 < z < 2.3. IChand et alJl2004 ISrianand et 811 12004). 

Strong constraints on variations in the electron-to- 
proton mass ratio, \x, ca n also be obtaine d from optical 
quasar absorption studies llThompsonlll975D . The first con- 
firmed molecular hydro gen quasar absorption line detec- 
tions were obtai ned bv | Lev shakoy fc Varshalovichl il985f) . 
and confirmed bv lFoltz et alJ Jl^8^L ^Iowever. although the 
large number of H2 transitions in the Lyman and Werner 



bands have different dependencies on /i, making them use- 
ful probes of /i-variation, unfortunately the relatively small 
fraction of DLAs in which H2 is detected has limited this 
work so far l|V arshalovich fc Leyshakov|l993|:|lvanchik et alJ 
l2002t lUbachs fc Reinholdll2004t llvanchik et alJl2005l) . Nev 
ertheless, the recent study of two high signal-to-noise spectra 
has yielded hints t hat n may have been smaller at redshifts 
between ~2. 5-3.1 jReinhold et al]l2006h \ 

A potentially more sensitive approach is possible when, 
as well as rest-frame UV, rest-frame 21-cm absorption from 
cold neutral hydrogen is also detected. In such cases, vari- 
ability of the constant x, defined above, can be investigated. 
Because the ratio of 21-cm to UV transition frequencies is 
proportional to x, it can be shown that, if both UV and 
21-cm absorption occurs at the same physical location, the 
relative change in the value of x between its value at redshift 
z and its value in a terrestrial lab is related to the observed 
absorption redshifts for rest-frame 21-cm and UV, Z21 and 
zuv, according to 



Ax/x • 



Xq 



ZlJV — Z21 



(2) 



X 1 + Z21 

(e.g. iTubbs fc Wolfelll980l) . In addition, because ^ enters in 
Equation [5] via the definition of x, results on the variation 
of /x may also be obtained, allowing completely independent 
tests of Afi/fi results obtained more directly from molec- 
ular absorption lines. However, application of this method 
is limited by the fact that there are only ~18 DLAs for 
which absorption has been observed both in the optical 
and in the radio llCurran et alJ 120051 . contains a detailed 
list). Consequen tly, there is only a handful of results based 
on this method JWolfe fc Davisill979l ITubbs fc Wolfdll980l 
ICowie fc Soneailal Il995t IWolfe et all Il98ll) . Among these, 
a high-resolution spectrum from a 10-m class telescope is 
used only once and for a single object dCowie fc Songailal 
Il995l) . Furthermore, there has been no attempt to apply 
this method to a sample consisting of more than one object. 
This is important for two reasons: 

(i) Any systematic effects due to spatial non-coincidence 
of hydrogen and heavy elements can only be detected if a 
sample of objects is used. 

(ii) To detect any space-time variation of x, one needs 
several observations, covering a range in redshift and direc- 
tion. 

The firs t results from su c h an approach were briefly 
presented in iTzanavaris et al.l (120051) . In this paper we ex- 
plain in detail how these results were obtained and extend 
our original sample of eight absorbers with a newly discov- 
ered 21-cm absorber. We have filled the gap in the litera- 
ture by applying this method to a sample comprising all the 
best available 21-cm absorption data in conjunction with the 
highest-resolution UV data available. The sample is made up 
of nine absorption systems situated along the line of sight to 
nine quasars covering the absorption redshift range ~ 0.23 
to ~ 2.35. The paper's structure is as follows: In Section |2] 
we describe our data sample. We describe our methods and 
results in Section [2] Our results are discussed in Section 0] 



1 Note that iReinhold et alJ i2006fl use the inverse definition of fi 
(i.e. fj, = m p /m e ) compared with this paper. 
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2 DATA 

2.1 Data sample and reduction 

Details of the data used in this work are given in Table Q 
and are discussed below. Our data sample is made up of 
nine quasar absorption systems along the lines of sight to 
nine quasars. Each system shows both 21-cm and UV ab- 
sorption (rest frame) along the line of sight to a different 
quasar. Rest-frame UV absorption is observed red-shifted 
in the optical region, and rest-frame 21-cm is observed red- 
shifted at longer radio wavelengths. In this paper we use 
the labels UV and 21-cm to distinguish the two types of 
observations that we use. 

Figures l3"lto II ll show the 21-cm and UV data registered 
to a common velocity scale. We note that all data have been 
corrected to the heliocentric frame. 



2.1.1 Radio data 

The original radio 21-cm absorption spectra were not 
available for this work, apart from the data for quasar 
Q0458— 020 which were kindly provided by Art Wolfe. 
Therefore, we digitised the 21-cm plots from the literature 
(see the references in T able using the dexter applica- 
tion llDemleitner et alJlioOlh . available at NASA/ADS for 
on-line scanned papers or, for PDF plots, the stand-alone 
version available at http://sourceforge.net/projects/dexter 

The initial step was to digitise each spectrum, obtain- 
ing flux/abscissa values. Then a linear wavelength scale was 
derived for, and assigned to, each spectrum by carefully mea- 
suring velocity or frequency points on the x-axis of the plot. 
This second step dominates the uncertainty in the zero-point 
of the velocity /frequency scale but is nevertheless not an 
important source of error. As an example, we digitised a 
postscript plot of the Sill 1808 A absorption feature ob- 
tained with Keck/HIRES and compared the resulting red- 
shift to that obtained from the original data. The two red- 
shifts differed by just 1 x 10~ 6 . A very conservative upper 
limit on the absolute velocity /frequency uncertainty is given 
by half of the digitised pixel size in each case. We checked 
the maximum contribution this zero-point error could make 
to the final estimate for Ax/x for each quasar and, in all 
cases, this was less than 25 per cent of the statistical-plus- 
systematic (i.e. increased) error bars illustrated in Fig. 
We therefore stress that the digitisation process generally 
contributes a small uncertainty in comparing the 21-cm and 
UV absorption redshifts. 

Redshift uncertainties and details of the observations 
and data reduction can be found in the references listed in 
Table However, since no redshift uncertainty is given for 
the rad io spectrum of Q1127— 145 inlKanekar^^ChenEalur 
i 20011 ) we used the uncertainty from lKanekar^^ClTBngahui 
i 2001aj) , as these are the same observers, using the same 
instrument to perform a similar type of observation. 

2.1.2 Optical (rest- frame UV) data 

Optical (rest-frame UV) high-resolution echelle spectral 
data for eight of the nine quasars were taken from the Eu- 
ropean Southern Observatory's (ESO's) archive. These were 
originally observed with the UVES spectrograph on the Very 



Large Telescope (VLT). D ata reduction was car ried out with 
the ESO/UVES pipeline .Ballester et aUl200(_) . modified to 
improve the flux extraction and the wavelength calibration 
precision. Post-pipeline processing was carried out using the 
custom- written code UVES popler 2 . For each quasar, the 
extracted echelle orders from all exposures were combined 
to form a single 1-dimensional spectrum, with cosmic ray 
signatures being rejected by a sigma-clipping algorithm dur- 
ing the combination. A continuum was constructed by iter- 
atively fitting a series of overlapping polynomials, reject- 
ing discrepant pixels at each iteration until convergence. 
For one of the quasars (Q1331+170) we have also included 
the part of the absorption spectrum for Sin 1808 A from a 
Keck/HIRES observation, provided by Art Wolfe. This spec- 
trum was reduced with the dedicated Keck/HIRES pipeline, 
makee 3 , written by Tom Barlow. 

For qu asar Q02 35+164 we digitised an absorption plot 
from lLanzetta fc Bowenl dl992T) for a single heavy element 
species. 



3 ANALYSIS 

There is one absorption system for each 21-cm/UV absorp- 
tion spectrum pair in a single quasar. To apply Equation 
we need values for zuv an< A z 2i- We now explain how we 
chose these values. 



3.1 Assumptions 

In each absorption system we used the frequency or wave- 
length of strongest absorption for each absorbing species, 
neutral or singly ionised, to provide the 21-cm and UV red- 
shifts. We justify these choices below but note that they 
need not represent the true physical picture for our results 
to be useful. Rather, they are used as a starting point for 
our analysis. 

3.1.1 Strongest absorption 

It is evident from Figs.l3T HT1 that in most absorption systems 
there is usually complex velocity structure. For neutral and 
singly ionised heavy element species there is sufficient sim- 
ilarity in this structure so that corresponding components 
can be identified. Unfortunately, this is not normally the 
case when one compares UV to 21-cm absorption profiles in 
the same system. In other words, it is not easy to tell which 
21-cm absorption component corresponds to which UV ab- 
sorption component. However, a decision must be made be- 
cause Equation [5] only holds if it is applied to 21-cm and 
UV redshifts for components which are situated at the same 
physical location. It is important to realise that this repre- 
sents a fundamental limitation in the comparison of 21-cm 
and UV absorption lines when deriving constraints on vary- 
ing constants. Clearly, this problem can only be overcome 
using a statistically large sample of absorption systems. The 
method used to determine the 'centroid' of the 21-cm and 

2 http:/ /www. ast.cam.ac.uk/~mim/UVES_popler. html 

3 http:/ /spider. ipac.caltech.edu/stafT/tab/makee/index. html 
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Table 1. Details of our data sample. There is one 21-cm/UV absorption system along each quasar sight-line. Column 1 is the quasar 
name and Column 2 its emission redshift. For the strongest component Column 3 gives the 21-cm absorption redshift. The uncertainty 
is given in parentheses and has been taken from the references shown in Column 4. The original 21-cm plots were also obtained from 
these references. However, the entries marked with f which refer only to the source of the 21-cm redshift uncertainty, whilst for the entry 
marked with * there is no uncertainty in redshift in the literature (see text). Column 5 gives the mean absorption redshift and standard 
deviation on the observed mean for the strongest UV component. Column 6 gives the UV heavy element species observed in the optical 
with the number of transitions in parentheses, if more than one. Column 7 gives the source for the UV data. For data obtained from the 
ESO/UVES archive, the ESO program identification code is given in this column and the principal investigator of the program is given 
in the footnotes. For Q1331+170 we also used Sill 1808 A Keck/HIRES data provided by A. Wolfe. For quasar Q0235+164 we digitised 
an absorption plot from the literature for a single heavy element species. 
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UV absorption is also clearly not crucial provided that it is 
unbiased and that uncertainties can be reliably estimated. 

Given the inherent dissimilarity between the 21-cm and 
UV absorption profiles, many possible methods exist for 
choosing appropriate 21-cm and UV redshifts. For example, 
for a single absorber, the various UV profiles might be cross- 
correlated with the 21-cm profile or one might attempt to 
identify 'nearest neighbour' UV and 21-cm velocity compo- 
nents. However, for the cross-correlation example, it is not 
clear how to derive a meaningful uncertainty on the redshift 
comparison. In the 'nearest neigh bour' approach, recently 
advocated in IKanekar et"al"l (2006), a detailed fit to the 21- 
cm and UV profiles is required to identify all the velocity 
components and, given the general variety of spectral reso- 
lutions apparent in Figs.lQ UTTl this is difficult. Moreover, this 
approach is strongly biased towards zero x-variation and this 
bias is again resolution-dependent. The redshift uncertainty 
for individual absorbers will therefore be too large. Simi- 
larly, for a large sample, the scatter in the values of Ax/x 
will be significantly underestimated and any additional scat- 
ter caused by intrinsic velocity shifts between the 21-cm and 
UV components will be masked. 

The simplest method is to take Z21 and ztjv from the 
strongest absorption component, as this is well-defined in 
our data. Although there is no general similarity in the ve- 
locity profiles, we assume that strongest absorption occurs 
at the same physical location for all species. We stress that 
this is a simple, initial assumption and later we demonstrate 
that we can actually use our results to test this assumption. 



Given the range of options above and given the fundamental 
limitation in associating 21-cm and UV components which 
are typically spaced by several km s , it is also clear that 
detailed fits to the components in order to decrease the cen- 
troid position uncertainty to the sub-pixel level is unimpor- 
tant. A cruder and simpler measurement, such as taking the 
redshift of the pixel with minimum intensity, will still involve 
an uncertainty which is smaller than the errors involved in 
comparing the 21-cm and UV profiles. 



3.1.2 Neutral and singly ionised species 

Neutral heavy element species might be considered more 
likely to be spatially coincident with 21-cm absorbing gas, 
which is due to neutral hydrogen. We decided to use singly 
ionised species as well for the following reasons: 

(i) In our sample, the velocity structure of neutral species 
was followed closely by singly io nised species. This wa s also 
true in the much larger sample in lMurphv et alj (2004) used 
to constrain variations in a. 

(ii) Extensive work by IProchaskal <l2003l) more gener- 
ally indicates that ionisation fraction, abundance and dust- 
depletion may not change significantly along an absorption 
complex. 

(iii) In any case, we only used the simple option of 
strongest components and these, in the overwhelming ma- 
jority of cases, corresponded very clearly between neutral 
and singly ionised species (see Figs. loT tTlTl . 
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By using both neutral and singly ionised species, we ob- 
tained a significantly larger sample of UV lines for each ab- 
sorber. This allowed us to investigate systematics for indi- 
vidual absorbers for the first time. Note though that there is 
one exception to our approach. This concerns Ca II which has 
not been used in our analysis, although it is singly ionised 
and shown in the figures for illustration purposes only (but 
see Section The reason for this is that the Can ionisa- 
tion potential, 11.871 eV, is such that one would not expect 
this species to trace neutral components well compared to 
the other singly ionised species used. Since DLAs have high 
neutral hydrogen column densities, one expects only photons 
with energies < 13.6 eV to penetrate the radiation shielding 
gas envelope. Such photons would then be unable to ionise 
species such as Fen (ionisation potential 16.18 eV) or Mnn 
(ionisation potential 15.64 eV), but they would still be able 
to ionise Can. Thus all singly ionised species used, apart 
from Call, are much more likely to be the dominant ionised 
species in an absorption system. For this reason we did not 
use this transition, just as we did not use higher ionisation 
species. 



3.2 Redshift determination 

In the spectral region covered by an absorption system, 
strongest absorption occurs at that point where the observed 
intensity is at a minimum. For each 21-cm absorption com- 
plex we thus identified and measured the dispersion coor- 
dinate, MHz or km s _1 , at the pixel of minimum intensity, 
from which we then obtained 221. We then searched the op- 
tical data for heavy element absorption features close to the 
redshifts where there is 21-cm absorption. We thus identified 
a number of UV absorption features, some of which were due 
to neutral species and some due to multiply ionised species, 
with the majority due to singly ionised species. For all de- 
tected UV features, where (1) absorption was due to either 
neutral or singly ionised species and (2) there was little or no 
saturation, we again determined the value of the dispersion 
coordinate, A or km s _1 , at the pixel of minimum intensity. 
We then determined the corresponding absorption redshift, 
zuv- When there were several UV transitions in a system, we 
carried out this procedure individually for each transition, 
e.g. independently for Znn 2026 and Znn 2062. In all, we 
used 31 different UV transitions. We present detailed veloc- 
ity plots of all 21-cm and UV absorption components used 
in this work in Figs. lo UTTl We include the Call profiles in 
these plots for comparison, even though we did not use them 
to determine zuv- These plots are centred at Z21, indicated 
by the solid vertical line. In each absorption system we cal- 
culated (zuv), the average of all redshift values determined 
for single species, as explained above. This is shown by the 
dotted vertical line in Figs. l3t flTl 



3.3 Ax/x Results 

We have obtained Ax/x results using three different meth- 
ods. The methods and results are summarised in Table [5] 
We have also performed a linear least squares fit to deter- 
mine the dependence of x on cosmic time. Additionally, we 
test evidence for any spatial variations across the sky. 



3.3.1 Method 1 

For each system we used (zuv) together with Z21 to ob- 
tain {Ax/x) aba, the average Ax/x value in a system, from 
Equation Since each quasar sight-line contains a single 
detected 21-cm absorber, each (Ax/x)^ value corresponds 
to a single quasar. These results are plotted in Fig.0 where 
labels are used to indicate which quasar corresponds to each 
plotted (Aa;/a;} a bs result. 

The average of all (Aa;/a;} a bs values is (Aa;/a;)totai = 
(0.52 ±0.99) x 10~ 5 . We label this 'result V. This covers the 
absorption redshift range 0.23 < z a b s ^ 2.35, correspond- 
ing to look-back times between 2.7 and 10.5 billion years 
4 . For this result we take as our quote d error the s tandard 
deviation on the mean result, (e.g. see lMevedll99l p. 23), 
i.e. 



n(n 



(3) 



where u = (Ax/a;)totai, iM = (Ax/x) ahs and n = 9, the 
number of quasars used. 

In other words, for method 1 we have not used individ- 
ual errors on Z21 and (zuv) but instead used the observed 
scatter to estimate the error on the mean. However, as ex- 
plained next, in methods 2 and 3 we have used the individual 
errors. It turns out that these further calculations provide 
an a posteriori justification of method 1 (see Section 



3.3.2 Method 2 

As mentioned, for Z21 we used the uncertainties given in 
the references listed in Tabled (see Section [2.1.11 . For zuv 
we used the standard deviation on (zuv) in each absorption 
system, i.e. 



E(w- 



n(n - 



(4) 



where w = (zuv), uli is a zuv measurement for a single 
transition and n is the number of UV transitions used for 
an absorption system. We note that it is not possible to de- 
termine o.^ in the cases of quasars Q0235+164, Q0952+179 
and Q0827+243, for which there is only one UV transi- 
tion (remembering that we do not use Call). Instead, in 
these cases we have used the value 6 x 10 -6 which repre- 
sents a typical error out of all the (zuv) values, except for 
Q1229— 021 and Q0438— 436 which are atypical cases. In the 
case of Q1229— 021, the larger error in (zuv) is due to the 
fact that the strongest component in Ti 11 3384 is located at 
~ +38 km s _1 from z%\. (see Fig. [KJ. This is the only case 
among the UV transitions used in all spectra where there is 
such a large individual discrepancy. However, for consistency 
and simplicity, we have used the redshift value obtained from 
the Till 3384 component to calculate (zuv). In the case of 
Q0438— 436, both the 21-cm spectrum and the UV spec- 
tra have quite low signal-to-noise ratio and dust-depletion 



4 We use a Hubble parameter Hq = 73 km s — 1 Mpc -1 , a total 
matter density !!m = 0.27 and a cosmological constant Qa = 0.73 
throughout this paper. The age of the Universe in this model is 
13.3 billion years. 
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-10 
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• This work 
□ Reinhold et al. (2006) 
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12 3 

Absorption redshift, z abs 

Figure 1. Ax/x results for the nine absorption systems in our 
quasar sample. Initial (increased) error bars have wider (nar- 
rower) error-bar terminators. Quasar names are given truncated 
to four digits. Each point represents {Ax/x)^ s obtained from 221 
and (*uv)i f° r au heavy element species in the UV quasar spec- 
trum, versus average (^tjv) f° r that spectrum. The long-dashed 
horizontal line is result 3 and the short-dashed lines show the ±ler 
range. Under the assumption of Aa/a = our results translate 
to constraints on Afi/fi and can be compa red with the direct con - 
straints from two H2-bearing absorbers in lReinhold et all l2f) 



effects may give different positions for the strongest compo- 
nent in different transitions. However, these cases and the 
associated larger uncertainties do not reflect the usual situ- 
ation for the sample as a whole. 

For each absorption system we therefore used (zuv), 
Z21 and their uncertainties to calculate the uncertainty in 
(Ax/x) a bs (clearly (Ax/x)^^ itself remains unchanged for 
each quasar) . This uncertainty is shown with the wider error- 
bar terminators in Fig. and was used to obtain an overall 
weighted mean result of {Ax/x)™^ tcd = (1.99T0.30) x 1(T 5 
('result 2'). 

3.3.3 Method 3 

For result 2, the x 2 per degree of freedom, z/, about the 



weighted mean is x l v 



well beyond what is expected 



based on the size of the individual error bars (i.e. \ 2 / v ~ !)■ 
Given the difficulty in associating 21-cm and UV velocity 
components, this provides motivation to quantify this addi- 
tional scatter by increasing the individual uncertainties on 
(A:r/:r} a bs by an amount s so as to force x 2 /v to unity. That 
is, the total uncertainty for each absorber is given by 



1 2 



+ S 



(5) 



A value of x 2 / u = 1 is achieved when s — 2.10 x 10~ 5 . 
The error bars now include the statistical uncertainties (as 
in result 2) and a systematic component which stems from 
the problem of comparing 21-cm and UV profiles. With 
these increased error bars, the weighted mean result is now 
<As/x)™ c t i f 1 Mod = (0.63T0.99) x 1(T 5 which we label 'result 
3\ The increased uncertainty for each absorption system is 
shown with narrower error-bar terminators in Fig. We 
argue below that this new error reflects much better the 



Table 2. Ax/x result summary. Column 1 gives the label used 
for the method and result. Column 2 is a brief description of 
the error determination method. Column 3 mentions the type of 
Ax/x total result. Column 4 gives the actual result. Result 3 is 
our fiducial constraint on Ax/x. 



Label 


Error type 


Ax/x type 


Result [10^ 5 ] 


1 


Standard error 


(Ax/x) total 


0.52 ±0.99 




on the mean 






2 


Formal statistical 




1.99 ±0.30 




propagation 




3 


Increased error-bars 




0.63 ±0.99 



inherent uncertainty in the determination of Ax/x and, as 
a by-product, it provides the first reliable estimate for the 
line-of-sight velocity differences between 21-cm and UV ab- 
sorption. We treat this is our main result. 



3.3.4 Variation in x with cosmic time 

Using the increased error bars above, we performed a linear 
least squares fit, (Ax/x) abs = — x/x (iib)abs, where (tib)abs is 
the average look-back time calculated from the (zuv) values 
in an absorption system. The best fit gives x/x — (—0.6 ± 
1.2) x 10~ lj yr _1 . Note that by choosing this fitting function, 
the intercept of which is zero, we implicitly assume that 
Ax/x = at z — 0. Although this may be corroborated 
by terrestrial experiments, it has not been checked in other 
locations corresponding to z = with the same technique 
used here (e.g. using 21-cm/UV absorbers elsewhere in our 
own galaxy). 



3.3.5 Angular variation in x across the sky 

If very large scale (~10Gpc) spatial variations exist, one 
might expect x to be different in different directions on the 
sky. We plot the distribution of Ax/x over Galactic coordi- 
nates (I, b) in Fig. |21 The grey-scale indicates (Ax/x)^ and 
the size of each point scales with the inverse of the error on 
(Ai/i) abs . 

There is no obvious dipole (or otherwise) angular vari- 
ation evident in Fig. [5] We used a x 2 minimisation algo- 
rithm to find the best-fit dipole in the angular distribution 
of Ax/x. A grid of directions in equatorial coordinates, (RA, 
DEC), was constructed and for each direction, a cosine fit 
to Ax/x as a function of the angular distance, 0, to each 
absorption system is performed. This yields an amplitude 
for the dipole, (Ax/x)d, defined by 



Ax/x(4>) = {Ax/x)ZlT cd + (Ax/x) d cos . 



(6) 



The best-fit dipole has an amplitude (Ax/x)a = (2.0±2.5) x 
10" 5 and is in the direction Pa^(RA, DEC) = (7.0hr,48°) 
or P A x(l, b) = (169°, 22°) as marked on Fig.0 The dipole is 
not significant, despite the fact that we identified the best-fit 
direction (see further discussion below). 
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360 



-3 2 4 

Ax/x = (A/Vm)* ^ [10- 5 ] 

Figure 2. Ax/x results for our nine absorbers plotted over 
the sky in Galactic coordinates. The grey- or colour-scale en- 
codes the (Ai/iJjIjs values shown in Fig.^with increased error- 
bars. The size of the data points scale with the inverse error-bar. 
The pole (Pa x ) an( i anti-pole (Aa x ) of the best-fit dipole vari- 
ation in x is marked and the amplitude of the dipole is given. 
It is not statistically significant. The CMB pole and anti-pole 
are marked for comparison. Assuming Act /a = 0, our results 
constrain Afi/u and can be compared with the more direct con- 
straints on Au/fi bv lReinhold et ah! l2006t) . Including these ad- 
ditional points, the pole moves to Pa/j = (190°, 35°) and the 
amplitude becomes (Au/fi)^ = (4.9 ± 1.7) X 10 -5 , which is for- 
mally significant. However, the dipole model is only preferred over 
the constant/monopole model at the 1.6-it level (see text). 



3.4 A/z//n results 

Results on Ax/x can be used directly to obtain limits on 
the variation of p. From the definition of x we obtain 



Mo 



Furthermore, 
constants is 



= A/i/M = Ax/x - 2Aa/a - Ag p / 



(7) 



9p 



on fundamental 



the dependence of 
fairly weak: Ag p /g p 

where Aqcd is the quantum chro- 



q ^ A(m,/A QCD ) 

(m,; /Aq C d) 



ijFlambaum et al.l l20 
modynamic scale. Therefore the g p term in Equation [7| can 
be neglected, so that there is a simple relation between 
A/i/pi, Ax/x and Aa/a. Note that relations between the 
variation of a, fundamental masses and Aqcd based on 



( 1984") : ICalmet & Fritzschl J2002I): iLaneacker et al. 


(■I 


Olive et ah] 




2002); 


Dent &i Fairbairnl 


(2003): iDine et alJ 



(2003) and references therein. Within these models the vari- 
ation of a is 1 — 2 orders of magnitude smaller than the 
variation of fundamental masses and Aqcd, so it may be 
neglected in the variation of x. However, these relations be- 
tween the variations of different fundamental constants are 
strongly model-dependent. 

Combining our result 3 for Ax/x with the claim ed de- 
tection of variations in a from iMurphv et alJ (120031) gives 
A/j,/h — (1.7 ± 1.0) x 10~ 5 . In this case, the non-zero 



Aa/a drives the resulti ng non-zero A / i/^t. U sing the claimed 
null-result on Aa/a of IChand et alJ ll2004l) yields Afi/fj, = 
(0.75 ± 1.0) x 10~ 5 . Although these two varying-a results 
are quite different, note that the error-bar in the resulting 
value of Afi/fi is completely dominated by the error in Ax/x. 
Alternatively, by conservatively assuming Aa/a — 0, we ob- 
tain A/1//1 = (Ax/x)Zi!t tcd = (0.63 ± 0.99) x TO" 5 . 

Note that the dependence of the proton mass on 
the current quark mass, m , is ve ry weak (Am p /rn p « 
0.05Am./m.. lFlambaum et alJliool) . To high accuracy the 
proton mass is proportional to the strong interaction param- 
eter Aqcd- Therefore, the limit on fi = m e /m p is simulta- 
neously a limit on the more fundamental Standard Model 
parameter v = m e /AQCD since Afi/ /x « Av/v. 



3.4-1 Variation in /i at high-redshift 

iReinhold et alJ ll2006l) used two HVbearing absorbers, one 
at z a bs = 2.594 towards Q0405— 443 and the other at z a b s = 
3.025 towards Q0347-383, to derive A/x/m = (-2.78 ±0.88) 
and (— 2.06±0.79) x 10 -5 , respectively, (formally) indicating 
a significant detection of cosmological variations in /1. Their 
quoted mean result is (—2.4 ±0.6) x 10 (weighted fit) and 
(-2.0±0.6) x 10" 5 (unweighted fit). Taken at face value, the 
error-bar on A/x//x from our results is comparable with the 
precision of these latest direct constraints on Afi/fi. In that 
case, even assuming that Aa/a — at high redshift, our 
null result 3 on (Ax / x)^^ ted appears to contrast sharply 
with these more direct measurements of Afi/fi. 

The Reinhold et al. data are plotted in Fig. Q for com- 
parison with our constraints on Ax/x. If the apparent dis- 
crepancy between these two data-sets is to be resolved by 
appealing to time-variations in \x, then clearly very sharp 
and possibly oscillatory variations would be required. How- 
ever, our full data set shows a significant gap in redshift, 
Az ~ 1.25. This naturally divides our set into a low- and a 
high-redshift subset. As the two Reinhold et al. datapoints 
are at high redshift, it is clear that they are best compared 
with our subset of four high-redshift points, the results for 
which are shown in Table As can be seen in this table, 
apart from the case of pure statistical propagation of errors, 
our results have a factor of ~ 2 lower precision than those 
of Reinhold et al., and, given the uncertainties, are, in fact, 
in agreement (2a). 

3.4-2 Variation in fi with cosmic time 

Using a similar approach to that for x/x above, we assume 
a linear variation of /1 with time to constrain \i/\i. With 
the increased error bars on the (Ax/x) & \, s values, assum- 



ing A a/a — and including the results of IReinhold et ail 
(2006), the time-variation of fi is constrained to be = 
(1.7 ±0.5) x lO^yr" 1 . Again, we assume that A/j,/ fi = 
at zero redshift but, once more, we emphasize that no 
absorption-line checks at z = have been made to justify 
this assumption. 



3.4-3 Angular variation in jx across the sky 

Figure [5] shows the two Reinhold et al. absorbers plotted 
on the sky with our A/1//1 values (i.e. our (Ax/x)^ values 
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Table 3. Afi/fi result summary. These results were produced 
in the same way as those tabulated in Table 1^1 but only the four 
highest-redshift absorbers have been u sed. Rows 1, 2 and 3 a are 
for Aa/a = 0, row 3b for Aa/a fro m|^hir phv et all 120031) and 
row 3c for Aa/a from lChand et alj i2004ft . 



Label 



Error type 



An/n type 



Result [10" 



1 


Standard error 
on the mean 


(A/V/i)total 


0.84 ± 2.00 


l 


2 


Formal statistical 


(^/t& cd 


2.04 ±0.31 


i 




propagation 


<A^>r l s i htcd 






3a 


Increased error-bars 


0.58 ± 1.95 


l 


3b 






1.73 ± 1.96 


2 


3c 






0.70 ± 1.95 


3 



l:Aa/a = ; 2: Aa /a from iMurphv et al.l <2003ft : 3: Aa/a from 
IChand et all i2004h 



assuming Aa/a — 0). Given the concentration of points with 
Afi/fi < below the Galactic equator and the mix of values 
above, it is tempting to explain the apparent discrepancy 
between these two data-sets by spatial variations in n across 
the sky. We performed a similar dipole search as for x above, 
including the Reinhold et al. data, finding the best-fit dipole 
direction to be Pa m (RA, DEC) = (8.5hr,33°) ar PA^(l,b) = 
(190°, 35°). The dipole amplitude is (Afj,/n) d = (4.9±1.7)x 
10 -5 . 

However, the limited range of angular separations in 
Fig. [5] severely limits the dipole interpretation. Although 
the dipole amplitude is formally significant at the 2.9-cr 
level, this is only after we had selected the best-fit direc- 
tion; this does not answer the question 'how significantly is 
the dipole model preferred over the constant (i.e. monopole) 
model?'. We address this question using a bootstrap tech- 
nique. Bootstrap samples are formed by randomising the 
values of Afj,/ ' [i over the quasar sight-line directions and the 
best-fitting dipole direction and amplitude is determined for 
each bootstrap sample. The resulting probability distribu- 
tion for (A/i/p)d indicates that values at >2.9-<j significance 
occur 12 per cent of the time by chance alone. Therefore, the 
data only marginally support significant angular variations 
in Afi/fj, given the sparse coverage of the sky. 



4 DISCUSSION 

The main result of this work is that the data show no evi- 
dence for any time variation of Ax/x. It can be seen from 
Fig. that the {Ax/x)^ points show significant scatter 
beyond that expected based on their statistical (i.e. origi- 
nal) error bars. The error in result 1 directly reflects this 
scatter. After taking into account individual redshift errors 
and propagating them statistically, we need to increase in- 
dividual uncertainties so as to obtain ax 2 /"~l about the 
weighted mean. Even so, this procedure (method 3) gives 
an overall error which is essentially the same as for result 
1. This simply says that the scatter in (Ax/x)^ completely 
dominates any errors in individual redshift measurements. 
This is an a posteriori justification of the validity of method 
1, in which individual redshift errors were neglected and not 
propagated statistically. The importance of this has not been 



realised before, as this is the first time that Ax/x has been 
determined for a sample of objects. If one determined Ax/x 
for a single quasar, there would be no scatter, but there 
would still be a hidden and unquantifiable systematic error. 
In practise, of course, for a sample of systems, as we have 
here, this additional error has random sign and magnitude 
and can be treated statistically. However, our result shows 
that any purely statistical determination of uncertainties, 
i.e. without taking the observed scatter into account, would 
inevitably underes timate the true unc e rtaint y. 

For example, ICowie fc Songailal dl995l) used neutral 
carbon lines (Cl and Cl*) observed in the Keck/HIRES 
absorption spectrum of Q1331+170, to obtain Ax/x = 
(0.70 ± 0.55 at at) x 10 -5 . For the same object we obtain 
Ax/x = (-2.59 ± 0.74 sta t ± 2.10 syst ) x 10~ 5 . The similar- 
ity in the statistical errors demonstrates our point. Here our 
statistical error is just 



(- 2 



< z >uv 



(1 + 221 



(8) 



and the systematic error is simply the value for the addi- 
tional uncertainty s which gives a % 2 l v value of 1 (Sec- 
tion l3.3.3p . We can only determi ne this because we have se v- 
eral objects. The central value of lCowie fe Songailal 1^)95) is 
different because they used a zyy value which is a wei ghted 
mean for observed components l|Songaila et alJ fl994ll and 
leads to a value very close to Z2i- Further explanation of 
the discrepancy in the central value may come from the fact 
that, in this system, we have used 23 distinct UV heavy 
element transitions with strongest components well defined 
within a few km s . 

The preponderance of the scatter over redshift uncer- 
tainties in determining the overall error also eliminates any 
worry about errors introduced via the digitisation process. 
We emphasise that having the original digital data would 
not have lead to appreciably different results in this work. 

What does this quantified scatter, s, mean? Simply 
from looking at the 21-cm and UV absorption profiles in 
Figs. IST tlll it clearly means that the assumption that 21- 
cm and UV absorption occurs at the same physical location 
is not correct. In other words the 21-cm and UV absorb- 
ing gases are randomly offset in space. A rough estimate 
of the average line-of-sight velocity difference is given by 
Allies ~ cs — 6.3 km s . We may explain this difference 
qualitatively. The emitting 21-cm quasar source may appear 
to have a large angular size to the absorber. A 21-cm sight- 
line can then intersect a cold, neutral hydrogen cloud with 
little or no heavy elements, whilst a UV/optical sight-line 
can intersect another cloud with heavy elements at quite a 
different velocity. Such a large angular size may be due to 
the combined effects both of physical proximity of the ab- 
sorber to the quasar and of intrinsic size of the radio emit- 
ting region. An alternative explanation involving small scale 
motion of the inte rstellar medium has been suggested by 
ICarilli et ail fcOOOl) in the case of the ~ 10 km s 1 velocity 
difference between 21-cm and mm absorption lines. 

The robustness of our Ax/x and Afi/fi results, as well as 
our estimate of Aui os , can be investigated by means of alter- 
native calculations. Table 2] summarises the simple tests we 
have carried out. All results quoted here have been obtained 
by applying method 3. Firstly, if we use 221 directly from the 
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Table 4. Alternative results. The results shown in Columns 2, 3 and 4 were obtained by introducing the changes briefly described in 
Column 1 (see text for more details). All results shown here were obtained using method 3. 



Change 


<^/<ol S i htCd A0- 5 


A/VAt/!0~ 5(a) 


A^/IO" 5 ^' 


Av loB 




quoted result 


0.63 ± 0.99 


1.78 ±1.02 


0.63 ±0.99 


6.3 km s~ 


-1 


221 from literature 


0.97 ± 1.05 


2.11 ± 1.07 


0.97 ± 1.05 


6.8 km s~ 


-1 


+ no digitised plots 


1.20 ±1.22 


2.35 ± 1.24 


1.20 ± 1.22 


8.0 km s" 


-1 


221 from Gaussians 


0.49 ± 0.94 


1.64 ±0.96 


0.49 ±0.94 


5.8 km s~ 


-1 


with Call 


0.71 ± 0.96 


2.38 ± 0.99 


0.71 ±0.96 


6.1 km s~ 


-1 




a Aa/a fromlMurDhv 


et al.N2003h: b A 


a/a = 0. 







litera ture, we obtain, for ( Ax / x)™^f^ tcd , A/1//1 (with Aa/a 
from iMurphv et all 120(531 and A/x//x (with Aa/a = 0), 
0.97 ± 1.05, 2.11 ± 1.07 and 0.97 ± 1.05, respectively, and 
in units of 10 (row 2 in Table 2J. If, additionally, we ex- 
clude any results for Q0235±164 and Q1127— 145, we obtain 
1.20±1.22, 2.35±1.24 and 1.20±1.22 (row 3 in TablegJ. We 
have also fitted Gaussians to the 21-cm absorption profiles 
and used these to calculate 221, i.e. we used the location of 
the minimum in the fitted Gaussian profile rather than the 
pixel of minimum intensity. We then obtained 0.49 ± 0.94, 
1.64±0.96 and 0.49±0.94 (row 4 in TablegJ. Finally, we have 
investigated what happens if one includes the Call lines in 
the calculations as well. This gives 0.71 ± 0.96, 2.38 ± 0.99 
and 0.71 ± 0.96. Note that, because we include Can, we 
have no problem calculating (Section 13. 3. 211 for quasars 
Q0235±164, Q0952±179 and Q0827+243. 



It is thus obvious that our results are relatively insen- 
sitive to the method used to estimate absorption redshifts, 
and, perhaps, even to the ionisation level of the transitions 
used. The dominant factor is the scatter due to the veloc- 
ity mismatch between 21-cm and UV transitions. Because 
the uncertainties in our determinations of UV redshifts are 
small when compared to this velocity mismatch, we do not 
expect that standard Voigt profile fitting of the UV absorp- 
tion features would significantly affect our results either. The 
mean line-of-sight velocity differences may become compa- 
rable in magnitude to individual errors in redshift measure- 
ments only if a much larger (~ 100) sample of 21 cm/UV 
absorbers is obtained: The average uncertainty per absorp- 
tion system due to the velocity offset is now ~ 10 but with 
100 absorbers, or so, we would expect this to be reduced by 
a factor of \/T00, i.e. to ~ 10~ 6 . A definitive answer to the 
question of whether Ax/x varies, or not, with cosmic time, 
will have to wait until then. 



Another important result is that, even when we con- 
servatively assume Aa/a = 0, our measurements of Ax/x 
imply constraints on Afi/ n which are independen t of the 
current possible detections of lReinhold et all l)2006l) . These 
authors find a combined (unweighted) value of A/j,/ /j, = 
(—2.0 ± 0.6) x 10 from two H2-bearing absorbers while 
we obtain a less direct constraint of Afi/fj, — (0.58± 1.95) x 
10" 5 . This new constraint, based on four high-redshift 21- 
cm/UV absorbers combined with the assumption of Aa/a 
= Ag p /g p = 0, serves as an internally robust, independent 
comparison to the H2 measurements. 
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Figure 3. Velocity plot for 21-cm and UV absorption towards quasar Q0952+179. The 21-cm data are in units of flux (Jy). The quasar 
continuum (1.4 Jy) has been subtracted. The UV data are continuum normalised. The solid vertical line at km s — 1 is at Z2i- The 
dotted vertical line is at (ztjv)- In this and subsequent plots Call is shown for illustration only but has not been used in the calculation 
of the plotted (zyy). 
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Figure 4. Velocity plot for 21-cm and UV absorption towards quasar Q1127— 145. The 21-cm data are in units of flux (Jy). The quasar 
continuum (6.3 Jy) has been subtracted. The UV data are continuum normalised. The solid vertical line at km s is at 221. The 
dotted vertical line is at (z\jy). 
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Figure 5. Velocity plot for 21-cm and UV absorption towards quasar Q1229— 021. The 21-cm data are in units of antenna temperature 
(K), where the system temperature was 45 K. The UV data are continuum normalised. The solid vertical line at km s -1 is at Z2i- The 
dotted vertical line is at (ztjv)- 
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Figure 6. Velocity plot for 21-cm and UV absorption towards quasar Q0235+164. The 21-cm data are in units of optical depth. The 
UV data are continuum normalised. The solid vertical line at km s _1 is at 221. The dotted vertical line is at (ztjv). 
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Figure 7. Velocity plot for 21-cm and UV absorption towards quasar Q0827+243. The 21-cm data are in units of flux (Jy) with the 
quasar continuum (0.9 Jy) subtracted. The UV data are continuum normalised. The solid vertical line at km s — 1 is at 221- The dotted 
vertical line is at (ztjv)- 
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Figure 8. Velocity plot for 21-cm and UV absorption towards quasar Q1331+170. Both the 21-cm and UV data are continuum 
normalised. The solid vertical line at km s" 1 is at z%\. The dotted vertical line is at (ztjv)- 
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Figure 9. Velocity plot for 21-cm and UV absorption towards quasar Q1157+014. The 21-cm data are in units of flux density (mjy) 
after subtraction of a polynomial baseline. The UV data are continuum normalised. The solid vertical line at km s _1 is at 221. The 
dotted vertical line is at (ztjv)- 
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Figure 10. Velocity plot for 21-cm and UV absorption towards quasar Q0458— 020. Both 21-cm and UV data are continuum normalised. 
The solid vertical line at km s — 1 is at 221- The dotted vertical line is at (ztjv)- 
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Figure 11. Velocity plot for 21-cm and UV absorption towards quasar Q0438— 436. The 21-cm data are in units of flux density (mjy) 
after subtraction of a polynomial baseline. The UV data are continuum normalised. The solid vertical line at km s _1 is at Z2i- The 
dotted vertical line is at (jztjv)- 



